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ABSTRACT 


Both  single  channel  and  array  signal  enhancement  tech¬ 
niques  have  been  applied  to  Rayleigh  waves  from  a  small 
Greenland  Sea  earthquake  recorded  at  LASA  and  13  LRSM  or 
Observatory  stations.  Analysis  of  individual  LASA  long 
period  recordings  indicated  that  the  matched  filter  increased 
S/N  by  more  than  6  db  over  the  mean  S/N  of  band  pass  filtered 
(15-50  sec  period)  seismograms.  Band  pass  filtering  increased 
the  mean  S/N  by  only  1.5  db  over  the  mean  for  unfiltered  seis¬ 
mograms  for  LASA.  Additional  improvement  in  S/N  from  beam¬ 
forming  both  band  passed  and  matched  filtered  traces  approached 
the  expected  increase  for  uncorrelated  noise  provided  the  in¬ 
tersensor  spacing  (mesh-size)  was  at  least  30  km.  Comparable 
improvement  was  obtained  for  both  single  channel  and  beam- 
formed  LRSH  data.  For  LASA  the  beam- formed  matched  filter 
S/N  for  13  stations  was  17  db  above  the  mean  S/N  for  the 
individual  band  pass  filtered  seismograms;  that  for  13  LRSM 
stations  was  15-16  db  above  the  mean  S/N  of  band  pass  filtered 
seismograms.  Beam- forming  matched  filter  seismograms  consist¬ 
ently  produced  S/N  values  7-9  db  above  the  S/N  for  beam-formed 
band  pass  filtered  seismograms.  The  effects  of  such  array 
parameters  as  number  of  sensors,  sensor  spacing,  and  aperture 
on  signal  enhancement  are  evaluated  for  this  event. 


INTRODUCTION 


The  objective  ot  this  study  is  to  compare  the  effect¬ 
iveness  of  simple  single  channel  and  array  methods  for  Rayleigh 
wave  enhancement  and  to  determine  the  dependence  of  these 
methods  on  array  parameters  such  as  number  of  sensors,  sensor 
spacing,  and  aperture. 

From  the  practical  viewpoint  it  is  important  to  know  not 
only  how  well  each  available  signal  enhancement  -  noise  reduct¬ 
ion  method  performs,  but  also  whether  or  not  it  can  be  used 
routinely  on  a  large  number  of  recordings.  Only  with  this 
information  is  it  possible  to  judge  the  "cost  effectiveness" 
of  a  particular  method  for  routine  application.  In  this  study 
we  have  attempted  to  evaluate  and  compare  the  simplest  approaches 
or  combinations  thereof,  for  surface  wave  S/N  (signal  to  noise 
ratio)  improvement.  These  results  then  form  a  basis  for  future 
evaluation  of  more  elaborate  surface  wave  enhancement  methods. 

The  most  complicated  approach  we  have  included  is  array 
summing  of  the  simple  least  squares  matched  filter  seismograms, 
each  of  which  is  generated  by  the  cross-correlation  of  a  Z  com¬ 
ponent  trace  with  a  known  signal  wave-form. 

In  addition  we  have  investigated  the  influence  of  certain 
array  parameters  on  the  effectiveness  of  each  method  with  a 
view  toward  eventually  determining  the  relative  merits  of  the 
LASA  long  period  array  versus  the  dispersed  LRSM  network  for 
observing  long  period  surface  waves  from  small  teleeeismic 
events.  The  results  presented  in  this  context  should  be  re¬ 
garded  as  quite  preliminary  since  they  are  based  on  the  analysis 
of  only  one  event. 

The  single  channel  methods  applied  to  both  LASA  and  LRSM 


and  (3)  a  combined  band  pass,  matched  filter.  The  array  methods 
applied  to  LASA  were  (1)  phased  sum  of  raw  seismograms,  (2) 
phased  sum  of  band  pass  filtered  seismograms,  (3)  phased  sum  of 
phase  equalized  seismograms,  (4)  phased  sum  of  matched  filter 
seismograms,  (5)  phased  sum  of  combined  band  pass,  matched 
filter  seismograms  and  (6)  matched  filter  of  phased  sum  seis¬ 
mogram.  For  the  LRSM  array,  only  array  methods  (2),  (4),  and 
(5),  above  were  applied.  In  the  matched  filter  analysis, 

Rayleigh  waves  from  a  large,  well-recorded  event  were  used  to 
search  for  the  surface  waves  from  a  smaller  event  from  the 
same  region . 

The  data  used  in  this  study  were  recorded  18  November 
I960,  at  the  Montana  LASA,  LRSM  stations,  and  other  observatories 
by  long  period,  vertical-component  seismometers.  The  LASA 
station  locations  are  shown  in  Figure  1.  Each  LASA  sensor  is 
located  at  a  subarray  center  position,  and  each  responds  to 
seismic  waves  of  various  periods  in  a  manner  similar  to  that 
shown  in  Figure  2.  Locations  for  the  LRSM  and  Observatory 
instruments  are  shown  in  Figure  3,  and  their  approximate  res¬ 
ponse  is  illustrated  in  Figure  4.  The  epicenter  and  distance- 
azimuth  data  for  the  Greenland  Sea  events  used  are  given  in 
Table  I . 

In  the  following  sections  we  outline  in  some  detail  the 
methods  of  analysis  and  procedures  applied,  present  the  salient 
results  obtained  for  this  event,  summarize  the  important  con¬ 
clusions,  and  recommend  several  items  deserving  further  invest i- 
gat ion . 


METHODS  OF  ANALYSIS 

Several  single-channel  and  multi-channel  (array)  methods 
were  employed.  All  are  characterized  by  their  simplicity, 
straightforward  practical  application,  and  reasonable  computer 
time  requirements.  In  this  section  each  method  is  presented, 
but  approaches  which  are  well-known  and  routine  will  not  be 
developed . 


Single  Channel  Methods 

Band  pass  filter.  Numerical  band  pass  filtering  is  so 
commonplace  that  no  development  is  needed  here.  In  this  study 
a  non-recursive,  phaseless  filter  (band  pass  15-50  seconds 
period)  was  applied  in  the  time  domain.  This  filter  response 
is  shown  in  Figure  5. 

Matched  filter.  The  development  of  the  matched  filter 
approach  in  the  time  domain,  and  the  justification  of  its 
application  to  the  detection  of  surface  waves  from  small  events 
were  given  by  Alexander  and  Rabenstine  (Reference  1).  Basi¬ 
cally  the  technique  amounts  to  searching  a  record  x(t)  for  a 
known  waveform  y(t).  It  is  assumed  that  x(t)  =  ay(t)  +  n(t) 
where  a.  is  constant  and  n(t)  is  a  random  process.  The  matched 
filter  output  is  essentially  a  cross-correlation  between  x(t) 
and  y (t ) . 


In  this  particular  application,  y(t)  is  the  surface  wave 
from  a  large  event  recorded  at  a  given  station.  The  time 
series  x(t)  is  a  seismogram  recorded  at  the  same  station  and 
containing  the  surface  wave  from  a  much  smaller  event  having 
as  nearly  as  possible  the  same  epicenter  as  the  large  event . 
The  matched  filter  takes  advantage  of  the  fact  that  the  two 
signals,  having  traveled  essentially  the  same  path  will  have 


experienced  the  same  phase  distortion  due  to  dispersion,  so 
that  cross-correlating  them  eliminates  the  propagation  effects 
on  phase,  regardless  of  how  complicated  they  may  be. 

The  reason  the  matched  filter  is  effective  in  enhancing 
dispersed  surface  waves  is  that  it  compresses  the  long,  dis¬ 
persed  wave-train  into  a  pulse  of  short  duration,  while  random 
noise  is  not  compressed.  Thus,  the  energy  density  of  the  signal 
on  the  matched  filter  output  compared  to  the  energy  density  of 
the  noise  is  increased;  this  amounts  to  an  increase  in  S/N. 

The  mathematical  arguments  which  show  this  are  presented  in 
the  Appendix. 

Whereas  the  techniques  derived  and  the  results  presented 
in  the  above  referenced  report  were  obtained  operating  in 
the  time  domain,  identical  results  can  be  obtained  at  a  con¬ 
siderable  savings  in  computer  time  by  operating  in  the  fre¬ 
quency  domain. 

The  matched  filter  output  is  given  by 


c  (t ) 
xy 


iU)t 

X(w)  Y*  (uj)e  duj 


where  X(uo)  ~  frequency  spectrum  of  the  test  seismogram 
Y*  (ui)  =  complex  conjugate  of  the  spectrum  of  the 
reference  signal 

Using  the  Cooley-Tuhey  method  for  fast  Fourier  transforms, 
implemented  by  McCowan  (Reference  3),  we  transform  y(t)  and 
x(t),  form  the  product  X(W)  Y*(u>),  and  inverse  transform  to 
get  the  matched  filter  seismogram.  In  addition  to  speeding 
up  the  standard  matched  filter  computations  this  approach  has 
the  advantage  that  (a)  differences  in  instrument  response  can 


be  accounted  for  and  any  band  pass  filter  or  other  shaping 
filters  can  be  applied  with  negligible  increase  in  computing 
time,  (b)  searching  for  a  total  of  L  reference  signals  in 
x(t)  requires  only  L+l  Fourier  transformations  if'  the  re¬ 
ference  spectra  are  saved,  and  (c)  the  spectra  are  available 
for  other  uses  such  as  measurement  of  radiation  pattern  as 
a  function  of  frequency,  Q  measurements,  phase  velocity  cal¬ 
culations,  and  excitation  spectra  as  a  function  of  magnitude. 

Revised  matched  filter  computer  programs  were  written 
to  implement  this  frequency  domain  approach. 

The  computation  time  is  proportional  to  (MLogM)  •  K  •  L 
where  M  =  sampling  rate  times  the  window  length  (M  <  4000  in 
the  present  program)  K-  number  of  channels,  L  *  number  of 
regions  or  reference  signals.  For  a  sampling  rate  of  1  pt/sec 
and  21  channels  the  frequency  domain  matched  filter  would  allow 
continuous  search  for  a  given  reference  waveform  at  about  6 
pts/sec.  This  means  that  one  could  search  21  seismograms  con¬ 
tinuously  for  events  from  approximately  six  different  source 
regions,  using  the  CDC  1604B  computer. 

Combined  band  pass,  matched  filter.  This  was  done  in 
either  of  two  equivalent  ways: 

a)  first  band  pass  filter  the  seismogram  and  then  use 

:his  filtered  seismogram  as  input  for  the  matched 
filter . 

b)  do  both  in  one  pass  through  the  frequency  domain 
version  of  the  matched  program. 

Multi-channel  Methods 


Phased  sum  of  raw  seismograms .  The  approach  used  in  beam¬ 
forming  is  so  well-known  that  no  elaboration  is  needed.  Standard 
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beam-forming  programs  were  used. 


Phased  sum  of  band  pass  filtered  seismograms.  The 
standard  beam  forming  programs  were  applied  to  the  band  pass 
filtered  seismograms. 

Phased  sum  of  phase  equalized  seismograms.  A  program  was 
written  to  determine  the  transfer  function  for  the  LASA  array, 
that  is,  the  phase  velocity  as  a  function  of  frequency.  The 
larger  reference  event  from  the  Greenland  Sea  (Table  1)  was 
used  for  this  purpose,  and  all  computations  were  carried  out 
in  the  frequency  domain.  (The  interpretation  of  the  dispersion 
in  terms  of  LASA  structure  is  being  carried  out  in  a  separate 
study,  but  the  observed  dispersion  is  all  that  is  needed  to 
phase  equalize  the  seismogram  to  one  element  in  the  array) . 

The  equalization  transforms  the  waveform  at  one  station  into 
an  expected  waveform  at  some  reference  station  in  the  LASA 
array.  The  frequency  domain  matched  filter  program  has  an 
alternative  option  to  perform  the  equalization  on  the  smaller 
event  using  the  observed  phase  velocity  with  frequency  from 
the  larger  event . 

Once  each  array  element  has  been  equalized  the  usual 
array  summing  is  applied. 

Phased  sum  of  matched  filter  seismograms.  The  standard 
beam-forming  programs  were  applied  to  the  matched  filter 
seismograms  using  the  apparent  velocity  of  the  beginning  of 
the  reference  signal  window;  this  assures  that  the  matched 
filter  signal  peaks  will  align. 

Phased  sum  of  combined  band  pass,  matched  filter  seis¬ 
mograms  . 

The  procedure  here  is  the  same  as  for  item  4  above, 
except  that  the  combined  band  pass,  matched  filter  seismograms 


are  beam-formed.  This  results  in  signal  enhancement  due  to 
the  combination  of  (1)  eliminating  noise  outside  the  signal 
band  by  band  pass  filtering  (2)  increasing  the  signal  energy 
density  relative  to  the  noise  energy  density  on  each  seis¬ 
mogram  by  the  matched  filter,  and  (3)  N  improvement  in 
signal  to  noise  by  beam- forming. 

Matched  filter  of  phased  sum  of  band  pass  filtered  seis¬ 
mograms  .  In  this  approach  the  individual  seismograms  are 
first  beam-formed  and  then  match-filtered,  using  as  a  refer¬ 
ence  signal  the  phased  sum  of  individual  reference  signals. 

It  can  be  shown  that  this  will  produce  less  signal  to  noise 
enhancement  than  the  phased  sum  of  matched  filter  outputs 
unless  the  seismograms  are  phase  equalized  to  a  common 
station  before  summing  in  which  case  the  results  should  be 
equivalent.  However,  over  an  array  which  does  not  strongly 
disperse  the  signals  relative  to  one  another,  this  method 
should  be  almost  as  good  as  phased  summing  matched  filter 
outputs.  It  ha9  the  advantage  that  it  takes  only  1/KI  (N  - 
number  of  sensors)  as  long  to  compute,  since  only  one  matched 
filter  output  must  be  obtained. 


The  question  arises  as  to  what  measure  of  improvement 
in  signal  to  noise  is  appropriate  for  dispersed  surface  waves. 
The  problem  is  different  from  that  for  body  waves  in  that  the 
surface  wave  energy  is  progressively  distributed  over  a  longer 
time  interval  with  increasing  epicentral  distance  because  of 
dispersion  during  propagation.  This  makes  time  domain  measures 
of  signal  to  noise  such  as  \  peak-to-peak  amplitude  divided 
by  the  RMS  of  the  noise  somewhat  misleading  and  unrealistic. 
However,  such  a  measure  is  acceptable  for  the  matched  filter 
output  which  is  not  dispersed. 

In  spite  of  the  danger  in  using  H  (pk  to  pk)/RMS  as  an 
absolute  measure  of  S/N,  values  of  this  S/N  estimate  before 
and  after  applying  each  method  do  provide  a  realistic  measure 
of  the  improvement  in  S/N  achieved.  Therefore,  we  elected 
to  adopt  this  as  the  standard  definition  of  signal  to  noise 
throughout  this  study,  namely 

S/N  b  \  (peak-tc-peak  amplitude)  (RMS  noise)"1  (2) 

Partly  this  was  done  for  convenience  in  that  existing  programs 
could  be  used,  and  partly  to  permit  eventually  the  direct  com¬ 
parison  of  surface  wave  enhancement  with  long  period  body 
wave  enhancement  using  array  processing. 

We  chose  the  peak-to-peak  amplitude  of  the  Rayleigh 
waves  by  searching  the  seismogram  for  the  extremum  values 
over  a  time  interval  approximately  corresponding  to  a  group 
velocity  window  from  4  km/sec  to  3  km/sec.  For  determining 
the  matched  filter  peak-to-peak  amplitude  the  search  was  made 
over  on  interval  of  100  seconds  centered  on  the  expected  signal 


arrival  time.  Clearly  this  method  breaks  down  when  the  signal 
level  falls  below  the  maximum  noise  peaks  (w  3  x  RMS).  For 
this  reason,  some  of  our  initial  S/N  estimates  may  be  too 
large  and  thus  the  mean  S/N  against  which  we  compare  may  be 
too  large.  This  makes  our  estimates  of  improvement  conser¬ 
vative. 

It  should  be  pointed  out  that  because  of  the  pulse-like 
shape  of  the  matched  filter  signal  (see  Figure  7  for  example), 
\  of  the  peak-to~peak  amplitude  was  consistently  about  3  db 
smaller  than  the  maximum  positive  signal  amplitude.  This 
suggests  that  for  the  matched  filter  S/N  estimates  it  would 
be  more  appropriate  to  use  the  maximum  signal  amplitude  in 
the  future,  or  equivalently,  increase  the  present  S/N  values 
by  about  3  db.  Thus,  our  estimates  of  S/N  enhancement  using 
the  matched  filter  are  always  conservative. 


In  all  cases  the  RMS  estimates  were  made  using  at  least 
a  30  minute  interval  immediately  preceeding  the  Rayleigh 
wave  arrival  (4  km/sec).  This  interval,  of  course,  includes 
the  body  wave  phases,  resulting  perhaps  in  a  slight  increase 
in  the  RMS  values  obtained.  This  means  that  (a)  our  RMS 
values  represent  an  upper  bound  on  the  true  RMS  of  the  winter 
noise  at  each  station,  (b)  all  S/N  estimates  calculated  are 
conservative,  and  (c)  comparisons  of  S/N  improvement  are  still 
valid. 


The  measure  of  S/N  improvement  we  adopted  is  given  by 
the  formula: 

S/N  enhancement  (db)  =  20  Log 


,  (S/N)  on  a  sum  trace) 
K 

i  l  (S/N)i 


(3) 


i-1 


-  9  - 
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K  =  number  of  elements  used  to  form  the  sum  trace 

th 

(S/N)^  =  signal  to  noise  ratio  on  the  i  band  pass 
filtered  seismogram 

Therefore,  we  are  trying  to  improve  on  the  mean  of  the  K 
individual  S/N  values  for  band  pass  filtered  data.  In  this 
study  we  found  that  the  mean  for  various  subsets  of  stations 
never  deviated  by  more  than  1  db  from  the  mean  S/N  for  all 
21  LASA  elements  (4.1).  These  values  are  given  in  Table  IV. 


DATA  PROCESSING 


The  flow  chart  shown  in  Figure  6  outlines  the  sequence 
of  steps  followed  in  processing  the  seismograms  to  obtain 
the  desired  outputs.  LASA  long  period  seismograms  are  re¬ 
corded  in  multiplexed  digital  form  at  800  bits-per-inch  with 
a  sampling  rate  of  5  points/second,  whereas  the  LRSM  and 
Observatory  data  are  recorded  in  analog  form  on  frequency- 
modulated  tapes.  In  all  cases  the  sampling  rate  of  LASA 
seismograms  has  been  reduced  to  1  point/second,  and  the 
analog  data  have  undergone  A/D  conversion  and  decimation  to 
this  same  sampling  rate.  A  standard  sampling  rate  is  re¬ 
quired  for  both  the  matched  filter  and  the  array  summing 
programs.  Our  choice  of  1  point/second  meets  this  require¬ 
ment,  keeps  computer  usage  reasonable  without  sacrificing 
time  resolution,  and  permits  frequency  analysis  for  all 
frequencies  below  . 5  cps . 

Surface  waves  from  the  events  were  sequentially  recorded 
at  the  LASA  on  two  multiplexed  tapes.  The  first  contained 
a  part  of  the  surface  wave  signature  of  the  earlier  event, 
and  the  second  held  not  only  the  remaining  data  from  the 
earlier  event,  but  also  the  entire  wave  train  from  the  later 
earthquake.  Sampling  rates  on  both  tapes  were  reduced  to  1 
point/second,  and  the  data  were  formated  and  merged  onto  a 
single  library  tape.  The  resulting  tape  was  used  to  form 
two  SUBSET  tapes,  one  containing  60  minutes  of  data  (Rayleigh 
minus  30  minutes  to  Rayleigh  plus  30  minutes),  and  the  other 
containing  10  minutes  of  the  Rayleigh  signature  from  the 
larger  earthquake  including  data  in  the  velocity  window  4  to 
3  km/sec.  At  this  point  we  deviated  from  normal  procedure 
to  remove  spikes  caused  by  the  tape  merger  mentioned  earlier. 


This  was  accomplished  by  a  special  program  which  vised  a  straight 
line  interpolation  method  to  replace  the  spikes.  As  a  result 
of  these  steps,  we  produced  two  LASA  SUBSET  tapes  which  were 
subsequently  used  as  input  to  the  matched  filter,  band  pass 
filter,  and  summing  programs. 

The  procedures  applied  to  the  LRSM  and  observatory  data 
were  essentially  the  same  as  those  applied  to  the  LAFA  data 
except  for  A/D  conversion  and  deletion  of  tape  merging. 


RESULTS 


In  this  section  we  present  the  major  results  obtained 
from  the  analysis  of  the  small  Greenland  Sea  event  documented 
in  Table  I.  All  noise  es  imates  given  represent  an  upper 
bound  on  the  RMS  levels  at  each  station  since  in  most  cases 
the  noise  sample  intervals  included  the  body  waves.  Values 
of  S/N  improvement  obtained  are  to  some  extent  systematically 
underestimated ,  (a)  because  estimated  signal  amplitudes  before 
processing  were  too  large  for  those  stations  where  the  maximum 
noise  peak  was  above  the  maximum  signal  level  and  (b)  because 
the  matched  filter  signal  estimate  we  use  is  always  less  than 
the  maximum  positive  matched  filter  signal  amplitude.  Thus 
all  the  enhancement  results  we  present  may  be  regarded  as  con¬ 
servative  in  assessing  the  usefulness  of  a  particular  approach. 

Comparison  of  Single  Channel.  Enhancement  Methods 

Band  pass  filter:  The  band  pass  filter  appears  to  have 
been  only'marginally  effective  in  reducing  the  overall  noise 
level  for  LASA  (see  Table  III  and  Figure  7) .  This  implies 
only  that  there  was  little  noise  outside  the  signal  frequency 
band.  The  mean  of  the  RMS  noise  levels  of  the  21  LASA  LPZ's 
was  reduced  only  1.5  db  by  band  pass  filtering  with  a  pass 
band  from  15  to  50  seconds  (period).  The  generalization  that 
this  will  always  be  the  case  should  not  be  made.  In  many  cases 
including  the  LRSM  seismograms  for  this  event,  the  noise  field 
can  be  much  stronger  in  4-7  sec  microseisms  and/or  above  50 
sec  microseisms,  in  which  case  band  pass  filtering  will  cor¬ 
respondingly  be  much  more  effective.  For  some  of  the  LRSM 
stations  band  pass  filtering  led  to  noise  reduction  of  as 
much  as  11  dh.  However,  the  mean  for  all  13  stations  was 


reduced  by  5.4  db. 


Matched  filter;  The  capability  of  the  matched  filter  to 
compress  the  surface  wave  to  a  pulse,  or  band  limited  impulse, 
as  described  in  the  Appendix,  proved  to  be  a  much  more  effect¬ 
ive  device  for  increasing  the  S/N  ratio  than  band  pass  filter¬ 
ing  (see  Figure  7  and  Table  III).  The  mean  S/N  of  the  individ¬ 
ual  matched  filter  outputs  for  21  LASA  LPS ' s  was  6.1  db  higher 
than  that  of  the  band  pass  filtered  surface  waves.  A  similar 
improvement  was  obtained  for  most  of  the  LRSM  stations  where 
the  signal  before  matched  filtering  could  be  accurately  measured 
abova  the  noise  (See  Figure  11  and  Table  III). 

Reliable  estimates  of  S/N  improvement  could  not  be  made 
for  those  channels  where  the  signal  was  weak  relative  to  the 
noise  because  the  method  for  determining  signal  amplitude  does 
not  work  when  the  signal  is  below  the  level  of  the  largest 
noise  peaks  (»  3  x  RMS).  However,  visual  examination  of  the 
matched  filter  output  for  these  cases  indicated  that  the  signal 
was  usually  enhanced  so  that  it  stood  out  above  the  noise  back¬ 
ground.  Thus,  significant  enhancement  must  have  occurred  in 
these  cases  as  well.  In  synthetic  cases  Alexander  and  Raben- 
stine  (Reference  1)  show  that  typical  teleseismic  surface  waves 

k 

can  be  detected  at  a  sing'jie  station  using  the  matched  filter 
to  levels  of  the  order  of  S/N  =  .35  Therefore,  for  these  low 
signal  levels  a  different  method  must  be  devised  for  estimat¬ 
ing  S/N  before  matched  filtering  if  reliable  estimates  of  S/N 
are  to  be  .obtained. 

When  the  S/N  values  before  filtering  are  in  error,  then, 
they  will  be  too  large,  so  the  estimates  of  S/il  improvements 
which  we  obtained  are  conservative. 
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Band  pass  and  matched  filter:  The  mean  S/N  of  the  21 
LASA  LPZ  seismograms  processed  through  both  filters  above  was 
increased  only  a  negligible  amount  (  <  0.1  db)  over  that  of 
the  matched-f iltered-only  seismograms  (see  Figure  7  and  Table 
III).  This  is  as  one  would  expect  since  the  matched  filter 
not  only  compresses  the  signal,  but  also  acts  as  a  band  pas3 
filter  having  a  response  given  by  the  amplitude  spectrum  of 
the  larger  event  (see  Appendix) .  For  the  LRSM  stations  the 
mean  S/N  using  both  filters  was  about  1  db  better  than  for 
the  matched  filter  alone  (see  Table  III). 

Typical  results  comparing  the  effectiveness  of  all  these 
single  channel  methods  for  LASA  are  shown  in  Figure  7.  They 
are  arranged  in  the  order  of  increasing  effectiveness  with 
the  matched  filter  results  clearly  the  best,  although  there 
is  little  difference  between  the  two  matched  filter  cases. 
Figure  11  shows  single  channel  matched  filter  results  for  3 
of  the  13  LRSM  stations;  the  top  pair  depicts  the  best,  the 
middle  3et  typical,  and  the  third  pair  the  poorest  of  the 
single  channel  results. 

Comparison  of  Multi-Channel  Enhancement  Methods 

A  comparison  of  the  different  methods  for  the  multi¬ 
channel  cases  is  complicated  by  the  fact  that  the  effect¬ 
iveness  of  each  technique  depends  on  a  number  of  parameters, 
such  as  sensor  spacing,  array  aperture,  number  of  sensors, 
combining  weights,  etc.,  in  addition  to  the  character  of  the 
signal  and  noise  on  each  channel.  We  shall  discuss  the 
various  techniques  making  comments  as  to  effects  of  some  of 
these  parameters. 

In  all  cases  Brennan  weighting  was  used,  i.e,,  (1)  each 
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seismogram  was  normalized  to  unit  noise  power  by  dividing  it 
by  the  RMS  of  a  selected  noise  window  within  it,  then  (2) 
each  seismogram  was  weighted  by  its  S/N  ratio.  Brennan 
weighting  emphasizes  the  better  channels  while  de-emphasizing 
the  noiser  ones  in  all  cases,  such  that  the  combined  S/N  is 
optimum  and  always  greater  than  or  equal  to  the  best  individual 
S/N  value. 

Phased  sum  of  raw  seismograms :  The  phased  sum  of  raw 
(unfiltered) ,  band  pass  filtered,  and  matched  filter  seis¬ 
mograms  for  LASA  each  gave  a  noise  reduction  approaching  the 
square  root  of  N  (number  of  sensors)  providing  that  the  inter¬ 
sensor  spacing  used  was  greater  than  30  km  (see  Figure  13)  and 
14).  Spacings  less  than  this  include  too  much  coherent  noise 
to  achieve  /ft”  reduction.  Spacings  greater  than  30  km  do 
not.  seem  to  produce  further  noise  reduction.  Hartenberger 
(personal  communication)  has  obtained  similar  results  for 
LASA  using  a  number  of  different  noise  intervals. 

It  is  to  be  expected  that  as  the  array  aperture  is  in¬ 
creased,  the  effects  of  dispersion  will  cause  some  signal 
loss  in  a  simple  summing  process.  However,  in  this  experi¬ 
ment,  using  apertures  up  to  the  full  diameter  of  LASA,  this 
loss  was  not  more  than  2-3  db  (see  Figure  15  and  16).  This 
is  probably  due  mostly  to  the  way  in  which  signal  was  measured. 
Since  the  signals  were  shifted  on  a  velocity  very  near  that 
corresponding  to  the  undispersed  Airy  phase,  there  was  little 
change  across  LASA  in  Airy  amplitudes,  which  was  essentially 
our  measure  of  signal  amplitude.  However,  there  was  almost 
certainly  some  cancellation  of  signal  energy  at  frequencies 
different  from  the  Airy  phase,  which  have  varying  velocities. 


This  can  be  seen  by  careful  comparison  of  trace  1  with  trace 
3  in  Figure  8. 

Phased  sum  of  band  pass  filtered  seismograms!  All  of 
the  above  comments  apply  to  this  process  as  well.  In  addition, 
it  was  found  that  over  LASA  the  S/N  of  these  sums  exceeded 
that  of  the  unfiltered  sums  by  more  than  H  db  in  only  one  case, 
that  being  the  array  of  5  sensors  with  10  km  spacing .  In  all 
cases  except  this  one,  the  noise  outside  the  filter  band  was 
incoherent  enough  to  be  cancelled  by  the  array  combining  just 
as  well  as  by  filtering.  In  the  special  case  mentioned,  the 
filtered  sum  was  1  db  better  in  S/N  than  the  unfiltered  sum. 
Figure  8  shows  the  results  for  21  LASA  channels  compared  to 
results  for  the  other  methods.  Figures  9  and  10  show  similar 
comparisons  for  N  =  5  and  9  respectively  for  inter-sensor 
spacings  of  10  and  20  km. 

The  phased  sums  for  different  values  of  N (number  of 
sensors)  are  shown  by  solid  dots  in  Figure  12  with  sensor 
spacing  (mesh  size)  as  a  parameter.  For  example  the  combinat¬ 
ion  A,  B,  (5  sensors)  and  A,  B,  C  (9  sensors)  correspond  with 
an  inter-sensor  spacing  of  10  km.  Table  II  gives  the  para¬ 
meter  approximations  for  the  sub-groups  designated  on  the 
graph.  While  all  the  values  fall  below  J~ N?  if  one  connects 
points  of  equal  mesh  size  (A,G  to  A,C,D  for  example)  the 
slope  is  approximately  the  same  as  for  J  N?’  this  suggests 
that  although  N  *  is  not  achieved  the  percentage  improve¬ 
ment  on  adding  additional  sensors  keeping  the  inter-sensor 
spacing  fixed  is  the  same  as  for  uncorrelated  noise. 


The  S/N  enhancements  on  beam-forming  band  pass  filtered 
seismograms  for  various  subsets  of  LASA  stations  in  addition 
to  all  21  are  given  in  Table  IV,  Figures  13  and  14  show  S/N 


gain  in  db  as  a  function  of  approximate  inter-sensor  spacing 
(mesh-size)  for  LAS A  using  fixed  valuta  of  N  (number  of  sensors) 
of  S  and  9  respectively.  In  each  instance  uniform  density  of 
stations  was  maintained  with  inter-sensor  spacing  defined  as 
the  average  mesh  size.  Also  nhown  in  these  Figures  is  noise 
reductionomv.Suramirej.lt  is  clear  that  a  mesh  size  of  at  least 
30  km  is  required  at  LASA  for  achieving  enhancement  (noise 
reduction)  near  that  expected  for  uncorrelated  noise. 

Sum  of  band  pass  filtered  phase-equalized  seismograms; 

The  results  with  phase  equalization  of  LASA  seismograms  were 
in  all  cases  within  db  of  those  obtained  with  simple  phased 
sums  of  band  pass  filtered  seismograms.  This  process  should 
be  expected  to  work  much  better  thar  simple  time  shifting  at 
large  apertures,  since  phase  equalization  eliminates  the 
effects  of  dispersion  across  the  array.  However,  aa  mentioned 
above,  this  dispersion  did  not  produce  a  si-;  j  if  signal 

change  (by  our  definition  of  signal  amplitudt.,  ,  apertures 

up  to  200  km,  and  thus  it  is  to  be  expected  that  the  phase 
equalization  would  not  work  much  better  than  simple  time 
shifting.  This  is  evident  from  a  comparison  of  trace  2  with 
trace  3  in  Figure  8. 

It  might  appear  that  phase  equalization  would  be  a  con¬ 
venient  technique  to  eliminate  the  effects  of  dispersion  across 
a  very  large  (continental)  array;  however,  there  are  several 
problems  involved.  Firstly,  one  needs  a  good  phase  velocity 
model  for  the  entire  region  involved.  Secondly,  because  of 
the  azimuthal  range  involved,  the  only  common  point  where  the 
wave  forms  should  match,  and  therefore  to  which  one  can  equal¬ 
ize  the  phase,  is  the  epicenter;  and  when  one  equalizes  to 
the  epicenter,  one  has  essentially  performed  the  matched  filter 
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operation. 

Phased  sum  of  matched  filter  outputs,  in  summing  the 
matched  filter  outputs  one  maintains  the  increase  in  signal 
level  achieved  in  the  matched  filtering  of  the  individual 
channels  and  adds  to  this  any  gains  due  to  noise  reduction  in 
array  combining . 

The  phased  sums  of  matched  filtered  seismograms  for 
different  values  of  N  (number  of  sensors)  are  shown  by  open 
circles  in  Figure  12  with  sensor  spacing  (mesh  size)  as  a 
parameter.  As  mentioned  previously  Table  II  gives  the  para¬ 
meter  approximations  for  the  sub-groups  designated  on  the 
graph.  In  each  instance  the  matched  filter  sums  are  7-9  db 
above  the  band  pass  filter  sums  for  a  given  N  and  sub-set 
of  sensors.  The  slope  on  connecting  points  of  equal  mesh 
size  (AD  to  ADE  for  example) is  approximately  the  same  as  for 
J 1?.  Just  as  for  the  band  pass  results,  this  suggests  that 
additional  sensors , keeping  the  spacing  fixed,  gives  the  same 
percentage  improvement  in  s/N  as  for  uncorrelated  noise. 

The  S/N  enhancements  on  be am- forming  matched  filter 
seismograms  are  given  in  Table  IV  for  various  sub-sets  of 
LASA  stations.  Figures  13  and  14  show  the  enhancement  achieved 
by  phase-summing  LASA  matched  filter  outputs  compared  to  hand 
pass  filter  phased  sums  for  sub-sets  of  5  and  9  sensor^  * 
spectively  as  a  function  of  the  average  inter-sensor' spacing 
(mesh-size).  In  Figure  13  it  is  clear  that  for  a  given 
spacing  the  matched  filter  phased  sum  is  7-8  db  above  that 
for  the  band  pass  filtered  phased  sum  (shown  in  the  same 
Figures),  with  the  8  db  values  at  the  larger  mesh  sizes.  For 
the  9  sensor  case  (Figure  14)  the  matched  filter  sums  are 
from  8-9  db  above  the  band  pass  filtered  sum,  again  with  the 
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larger  differences  at  the  larger  spacings.  Since  the  noise 
reduction  values  are  very  similar  for  both  phased  sums,  the 
systematic  increase  in  the  difference  with  spacing  and  number 
of  sensors  probably  is  due  to  the  fact  that  because  of  dis¬ 
persion  across  LASA  some  signal  is  lost  in  phase— summing  band 
pass  filtered  seismograms  for  more  widely  separated  stations, 
while  little  is  lost  in  phase  summing  matched  filter  seis¬ 
mograms.  That  is,  aperture  affects  the  phased-sum  of  band 
pass  filtered  seismograms  more  than  it  does  the  matched  filtered 
sums.  This  is  shown  in  Figures  15  and  16.  Aperture  over  LASA 
has  little  effect  on  tiie  combination  of  matched  filter  outputs 
since  their  undispersed,  autocorrelation-like  wave  shapes 
coincide  exactly  as  long  as  they  are  shifted  properly.  Only 
significant  azimuthal  differences  in  the  source  phase  spectra 
between  the  reference  event  and  the  test  event  will  result 
in  matched  filter  signal  cancellation  on  summing  (see  Appendix) , 
whereas  with  the  previous  array  sums,  any  azimuthal  variation 
in  the  source  phase  spectrum  will  produce  some  signal  cancel¬ 
lation  on  summing  in  addition  to  that  due  to  dispersion  across 
the  array. 

One  might  expect  that  errors  in  the  epicenter  location 
or  origin  time  would  cause  difficulty  in  shifting  the  matched 
filter  outputs.  It  can  be  shown,  however,  that  if  the  seis¬ 
mograms  from  both  events  are  time  shifted  on  the  same  velocity, 
these  errors  tend  to  cancel  out.  In  fact  the  positions  of  the 
maxima  in  the  matched  filter  outputs  can  be  used  to  estimate 
these  errors. 

The  main  problem  involved  in  using  the  matched  filter  is 
the  compilation  of  a  library  of  master  events  to  cover  the 
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epicentral  areas  of  interest. 

Phased  sum  of  band  pass  and  matched  filter  seismograms. 

As  was  the  result  with  the  single  channel  seismograms,  sums 
of  LASA  seismograms  passed  through  both  filters  were  iden¬ 
tical  to  those  of  the  matched-filter-only  sums,  within  0.5 
db.  Thus,  all  comments  about  matched  filter  sums  are  appli¬ 
cable  here;  and  Figures  12,13,14  also  represent  this  method 
for  LASA;  It  can  be  concluded  that  in  most  cases,  band  pass 
filtering  before  matched  filtering  is  unnecessary  for  LASA. 

For  the  LRSM  stations  the  phased  sum  of  the  band  pass  matched 
filter  seismograms  was  better  by  about  1.4  db  than  for  matched- 
filter-only  sums,  so  that  band  pass  filtering  the  data  would 
seem  to  be  desirable  for  LRSM  stations  before  matched  filter¬ 
ing. 

This  method  gave  the  best  results  of  all  the  approaches 
we  used.  S/N  values  on  beam-forming  were  consistently  7-9 
db  above  the  beam-formed  band  pass  filtered  S/N  values  as 
can  be  seen  in  Table  IV.  LRSM  phased  sums  of  matched  filtered 
seismograms  gave  enhancements  about  the  same  as  LASA.  For 
13  sv,  isors  S/iJ  improvements  over  the  mean  of  individual  S/n 
values  was  17  db  for  LASA  compared  to  a  value  of  17.5  db 
expected  for  uncorrelated  noise;  for  13  LRSM  stations  the 
improvement  was  15  (16)  db  over  the  mean  of  individual  S/N 
values  compared  to  a  value  of  15.5  (16.5)  db  expected  for 
uncorrelated  noise. 

Matched  filter  of  phased  sum  of  band  pass  filtered  seis¬ 
mograms  .  This  was  tried  only  for  the  LASA  9  sensor  case  with 
30  km  spacing;  the  resulting  S/to  was  about  1-2  db  lower  than 
phase  summing  individual  matched  filter  seismograms.  This 
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suggests  that  if  this  much  loss  in  s/N  by  summing  before 
matched  filtering  can  be  tolerated,  this  approach  might  be 
preferable  because  of  the  savings  in  computing  time;  this 
method  takes  only  about  ^  (N  e  number  of  sensors)  as  long 
as  match  filtering  individual  channels  and  then  phase  sum¬ 
ming*  However,  the  reason  this  method  worked  as  well  as  it 
did  is  because  all  of  the  signals  were  similar  in  waveform* 

As  aperture  increases  and  the  dispersed  signals  become  dif¬ 
ferent  in  waveform  this  method  should  become  progressively 
more  inferior  to  the  previous  method.  Also  the  weighting 
becomes  a  problem  in  applying  this  approach,  since  we  do  not 
have  the  individual  S/N  values  for  the  matched  filter  seis¬ 
mograms  ;  the  Brennan  weighting  we  used  in  the  previous 
method  can  be  shown  to  produce  the  optimum  S/N  enhancement 
if  the  noise  is  uncorrelated  and  stationary.  Therefore,  the 
combination  of  dissimilar  waveforms  and  incorrect  weighting 
may  make  this  method  less  desirable.  At  this  point,  however, 
we  have  too  little  data  on  this  approach  to  assess  its  merits 
fully. 

Our  results,  then,  indicate  that  of  the  methods  investi¬ 
gated,  phase-summing  individual  matched  filter  outputs  pro¬ 
duces  the  greatest  enhancement  in  S/N  and  that  in  all  types 
of  phase  summing  of  long  period  data,  sensor  spacings  of  at 
least  30  km  are  needed  to  assure  that  the  noise  in  the  15-50 
second  period  range  is  uncorrelated. 


CONCLUSIONS 


For  a  single  Greenland  Sea  earthquake: 

In  By  means  of  a  matched  filter,  the  mean  signal/  noise 
ratio  for  the  surface  wave  on  21  LASA  LPZ  seismograms  was  in¬ 
creased  6  db  over  that  of  the  band  pass  filtered  seismograms. 
Mean  S/N  improvement  of  3  to  4  db  was  obtained  for  13  LRSM 
stations . 

2.  The  signal/noise  ratio  of  the  matched  filtered  seis¬ 
mograms  was  independent  of  whether  the  seismograms  were  pre¬ 
filtered  with  a  band  pass  filter  for  LASA.  Pre-filtering 
LRSM  seismograms  produced  matched  filter  results  about  1.5  db 
better  than  not  filtering. 

3.  An  additional  increase  of  signal/noise  approaching 
j  N  (N  *  number  of  sensors)  was  achieved  by  phased  summing 
the  matched  filter  outputs  for  LASA  if  an  inter-sensor  spacing 
of  at  least  30  km  was  maintained.  A  similar  improvement  was 
observed  for  the  LRSM  stations  which  had  a  still  larger  (but 
not  uniform)  spacing. 

4.  For  array  apertures  as  great  as  the  full  diameter 
of  LASA,  phased  equalized  summations  showed  little  increase 

(  <  1  db)  in  signal/noiae  over  simple  phased  sums,  both  having 
been  band  pass  filtered. 

5.  Phased  sums  of  matched  filter  outputs  were  consistently 
7-9  db  above  corresponding  phased  sums  of  band  pass  filtered 
seismograms . 

6.  A  comparison  of  matched  filter  phased  sums  for  13 
LASA  and  13  LRSM  stations  (spacing  >  30  km)  showed  S/N  gains 
of  17  and  15  (16)  db  respectively,  over  the  mean  of  individual 
band  pass  filtered  S/N  values.  In  both  cases  this  was  within 
h  db  of  the  value  expected  for  uncorrelated  noise. 


7.  Aperture  at  LASA  causes  little  or  no  signal  1-oss-  for 
matched  filter  phased  sums  and  only  moderate  signal  loss  (.5  to 

3  db)  on  hand  pass  filter  phased  sums  for  apertures  up  to  200  km. 
There  was  also  little  or  no  signal  loss  on  phase  summing  the  LRSM 
matched  filter  seismograms  over  a  continental  size  aperture. 

8.  Even  for  the  sensor  spacings  at  which  S/N  gains  were 
below  those  expected  for  uncorrelated  noise,  the  percentage  in¬ 
crease  in  S/N  on  adding  additional  sensors  was  approximately  the 
same  as  for  the  uncorrelated  case. 
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Several  aspects  of  this  study  require  further  work  and 
there  are  some  follow-on  areas  of  investigation  which  seem 
worthwhile.  Specifically  we  recommend: 

1.  Additional  small  events  be  studied  along  the  same 
lines  as  in  the  present  study  to  confirm  the  findings  of  this 
study  with  regard  to  S/N  gains  and  array  parameters. 

2.  An  attempt  be  made  to  define  experimentally  threshold 
signal  detection  levels  for  the  matched  filter  method  to  verify 
the  threshold  valut?s  previously  obtained  in  synthetic  cases. 

To  do  this  a  new  method  for  measuring  initial  S/N  values  must 
be  developed. 

3.  Further  work  be  done  to  assess  the  relative  merits  of 
matched  filtering  phased  sum  seismograms  versus  phase-summing 
individual  matched  filter  seismograms. 

4.  Attempt  to  use  both  single  channel  and  multi-channel 
matched  filter  approaches  to  obtain  AR  (ER)  estimates  for  weak 
signals. 

5.  The  empirical  determination  of  scaling  functions  to 
adjust  the  matchad  filter  reference  spectrum  to  that  appropriate 
for  the  event  of  interest. 

6.  Assessment  of  the  "cost  effectiveness"  of  the  maximum 
likelihood  matched  filter  approach  (Reference  1)  for  LRSM  and 
LASA  surface  waves.  This  will  at  the  same  time  give  additional 
information  on  the  characteristics  of  long  period  noise  at  the 
various  stations. 

7.  Further  experimental  work  along  the  lines  of  Reference 
1  and  the  Appendix  of  this  report  to  obtain  source  parameters 
such  as  radiation  pattern,  initial  phase,  and  excitation  spectra 
for  surface  waves  from  small  events. 


8.  An  attempt  be  made  to  detect  long  period  body  waves 
by  a  matched  filter  approach  where  the  reference  signal  consists 
of  several  body  wave  arrivals  spanning  a  time  interval  long 
enough  so  that  the  noise  in  that  interval  is  random  and  station- 
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APPENDIX 


We  present  here  a  discussion  of  the  details  of  the  matched 
filter  operation  which  shows  that  because  the  signal  is  com¬ 
pressed  in  time  while  the  noise  is  not,  S/N  is  increased „ 

The  test  seismogram  x(t)  has  the  Fourier  representation 


x(t)  -  /  X  (  to  )  e^^dto 


(a-  :l  ) 


and  energy  (by  Parseval's  identity) 


E  =  /  |  X(io)  dw 
—  00 


(A- 2) 


The  matched  filter  has  the  Fourier  representation 


c (t)  =  S  X (to)  y*(u)  eiwtdw 


(A- 3) 


where  Y*(w)  is  the  complex  conjugate  of  the  Fourier  transform 
of  the  reference  signal  y  (t) ,  (Y)  =  |  Y  (  w  )  I  )•  The  asso¬ 

ciated  energy  of  the  matched  filter  output  is 


Ec  =  /  |  X  (co) 
—  00 


I  Y  (to ) 


dio 


(A- 4) 


'J’his  ic  equivalent  to  the  energy  in  x(t)  after  it  is 
filtered  with  a  function  whose  spectrum  is  Y  ( to  )  .  Note  that 
if  v.e  whiten  Y(  (o )  (i.e.  divide  by  |  Y  (  <o  )|  )  so  that  instead  of 
Y  ■  (  w )  we  have  Y*  /  I Y (  w ) |  =  e  ^  in  equation  (A-3),  this 

equation  becomes 


cn(t)  -  Jx(»>e  -i[e(u,)  ‘  "t:!  dio  (A- 5) 

with  energy 


E  *  f  I  X  (u>)  1 2  dm  (A-6) 

cn  J  1  1 

_*> 

Thus,  if  the  reference  signal  is  whitened  before  apply¬ 
ing  it  as  a  matched  filter  the  energy  of  x(t)  is  conserved 

in  the  operation  since  E  *  E  . 

cn  x 

What  is  left  to  show  is  that  the  signal  is  compressed 
in  time  while  the  noise  is  not. 

If  we  talce  x(t)  =  s(t)  +  n(t)  (signal  +  random  noise), 

then , 


X(w)  -  I S (tu)  J  exp[i^cpos(w)  -  Jfa-  +  cp.  (■,)] 

+  j N ( u) )  |  exp^Y(m)  +  cpi(u>)^)j  (A_?) 

Y(oi)  -  |Y(w)  I  exp£-i^cpoy(u))  -  -~y  +  cp,(u)Q] 

■  |x(m|  e-ie(u,)  lA’8) 

and 

where 

cp  (<«)  =  phase  spectrum  at  the  source  for  s(t) 

O  3 

<Poy(tu)  ■  phase  spectrum  at  the  source  for  y(t) 

<Pt (a»)  *  instrument  phase  spectrum  (assumed  identical  for 

x  anu  y) 
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C  ( us )  =  phase  velocity  spectrum 

A  *  epicentral  distance 

Is (iu)  I  *  amplitude  spectrum  of  s(t) 

| Y(W) |  =  amplitude  spectrum  of  y(t) 

9  (to)  =  total  phase  associated  with  y(t) 
y(to)  =  random  phase  associated  with  n(t) 
Therefore, 


XY*  =  |S(U))| 

|Y(0))  j 

exp[i^cpos(ui)  -  cpoy(uj) 

N 

A 

+  In  (hi)  | 

|y(uj)  | 

expj^i^Y  (to)  -  6  (to)  ^ 

(A-9) 

The  first  term  on  the  right  hand  side  of  the  equation  (A-9) 
is  the  matched  filter  signal  3pectrum  and  the  second  term  is 
the  noise  spectrum.  We  v  ?.V;  discuss  each  term  separately.  From 
the  first  term,  the  signal  on  the  matched  filter  seismogram  is 
given  by 

00 

sf(t)  =  J  |  S  (tu )  |  |  Y ( ui )  |  exp£i{jPos(u>)  -  Voy(w)  +  (Utjj  dm  (A-10) 
—  0 0 

If  both  s(t)  and  y(t)  have  the  same  phase  spectrum  at  the 

source  then  =  «P  and  (A-10)  reduces  to 

os  oy 

sf  (t)  “  J  JS(UJ)  |  [  Y(uj)  J  eiUJtduj  (A-ll) 

00 

or  for  the  whitened  reference  signal  case 
00 

Sfn(t)  =  J  |a<<«)|e1U)tdu)  (A-12 ) 

_0Q 


A- 3 


t-  'l 


"P| 


-s  -i&rt'-.r;- 
f-*-. 


>-5.-  .-•»'  •"•  : 


/T  ■,'v 


We  are  now  in  a  position  to  compare  the  energy  density  of  a^(t) 
with  that  of 


s(t)  *  J  j  S (to )  |  expj^i  (jPoa(u>)  -  —  +  cpi(m)  +  uitjj  dU)  (A-13) 

—  00 

For  a£n(i)  the  maximum  energy  density  occurs  at  t  *  0,  since 
it  is  here  that  all  the  frequencies  add  together  in  phase  (zero 
group  delay  for  all  frequencies) and 


jjs(u))|dou  >  £  |  S  (tu)  |  eiU)tdu)  for  all  t  f  0. 


For  s(t)  the  energy  is  spread  over  a  longer  time  window  because 
the  group  delays  (energy  arrival  times)  are  frequency-dependent? 
that  is,  the  energy  arrival  at  frequency  U)  occurs  at  a  value  of 
t  such  that 


air  [<poa(ui)  -  ftl7  +  V®’  +  mtj  ■  »;a<u,)  -  UT?T  +  <pi,tu)  +  *•  ■  0 


A 

U(u») 


-  q>'  (w)  -  cp :  (uj ) 

OS  1 


(A-14) 


where 


U(w)  “  group  velocity  with  frequency 

cp1  (ou)  =  group  delay  time  at  the  source 
os 

cp^(tu)  *»  group  delay  time  through  the  instrument 

Since  this  total  group  delay  time  (A-14)  is  frequency 
dependent,  energy  at  different  frequencies  arrives  at  different 
times  on  the  seismogram.  Given  that  the  total  energy  of  s(t) 
is  equal  to  that  of  sfn(t)  by  the  same  arguments  presented 
earlier,  this  means  that  the  maximum  energy  density  of  s(t) 


.  *  .  .4.  Tin  , 


must  always  be  leBS  than  that  for  s^n(t)  (the  case  where  the 
energy  for  every  frequency  arrives  at  the  same  time).  Thus, 
the  signal  is  compressed  in  time  without  lose  of  energy-. 

We  now  show  how  the  noise  is  affected  by  the  matched 
filter  operation.  What  we  compare  is  the  noise 


00 

n(t)  *  j  |n(ui)  |  exp^ifjy(ui)  +  cp^(u))  +  u>tjj  <30U 


(A-15) 


with  the  matched  filtered  noise 


nf(t)  -  j  N(uu)  Y*(u))e1U)t  dtu 
—  00 

or  for  the  whitened  reference  signal  case 


(A-16) 


uu  A 


(t )  58  j  J N (u) )  |  exp^i^«PQy  (tu)  +  +  *(«»)  +  Wt^j  du)  (A-17) 


The  group  delay  at  frequency  UJ  in  (A-17 )  is  given  by 


t  (u»)  =  -r  (IB)  - 


uTUT  +  voyM 


(i-18) 


Since  ¥(uj)  is  random,  its  derivatives  ¥'(«)  are  also  random. 

The  second  two  terms  on  the  left  side  of  (A-18)  introduce  a 

systematic  shift  in  group  delay  time  with  frequency,  but  because 

¥ 1  ( uu )  is  random,  t^fui)  is  still  random  with  the  result  that  the 

noise  is  not  compressed.  Since  the  matched  filter  operating  on 

n(t )  in  the  whitened,  reference  signal  case  conserves  energy 

(by  the  earlier  arguments)  and  is  a  linear  filter  with  a  flat 

amplitude  spectrum,  the  noise  has  the  same  mean  (zero)  and  the 

same  variance  (RMS  amplitude)  as  before,  so  that  the  mean  energy 

density  in  n^  (t)  is  the  same  as  the  mean  energy  density  in  n(t) . 
f  n 


(See  Papoulis,  1965,  pp. 345-347  for  a  proof  of  this  point. 
Stationarity  is  assumed) . 

Thus,  we  conclude  that  the  matched  filter  operation 
increases  the  signal  energy  density  while  not  increasing 
the  noise  energy  density.  This  means  an  enhancement  in  any 
conventional  time  domain  S/N  estimate. 

Since  the  case  where  the  reference  signal  spectrum  of 
the  matched  filter  is  not  whitened  (the  usual  case)  is  equi¬ 
valent  to  prefiltering  x(t)  with  a  phaseless  filter  whose 
spectrum  is  ]Y  (u))  |  and  then  using  the  whitened  reference 
signal,  all  the  arguments  above  with  respect  to  increasing 
signal  energy  density  relative  to  the  noise  energy  density 
hold  for  this  case  as  well. 

Since  |  S(u>)  |  is  not  identical  to  |Y(U))|in  practice,  one 
should  shape  | Y ( to )  |  so  that  it  duplicates  |s(uj)|  as  nearly  as 
possible.  This  filter  |z(ui)i  | Y (tu )  j  ■  |S(u>)|  has  the  effect 
of  emphasizing  the  frequencies  where  signal  is  present  and 
de-emphasizing  the  frequency  ranges  where  signal  is  low.  As 
yet,  we  do  not  have  empirical  shaping  functions  |z(tu)|  as  a 
function  of  body  wave  magnitude,  so  for  the  present  we  simply 
use  |  Y  ( to )  |  .  Apparently  I Y ( to )  |  is  a  reasonable  first  approxi¬ 
mation  for  |  S(tu)  j  ,  because  the  matched  filter  results  for 
this  case  are  better  than  for  the  whitened  case  in  the  tests 
we  have  made  thus  far. 

As  a  final  note  it  should  be  pointed  out  that  all  the 
analysis  done  in  this  Appendix  assumes  infinite  limits  in  both 
time  and  frequency,  while  in  practice  we  work  with  finite  times 
and  finite  frequency  bands.  However,  it  can  easily  be  shown  that 
the  results  of  this  section  apply  equally  well  to  the  actual 
situations  involving  finite  intervals. 
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Table  III-A. 

Observed  S/N  Values  for  LASA  Stations 


STATION 

S/N 

Raw 

Seismograms 

S/N 

Band  Pass 

Filtered 

Seisroograirm 

S/N 

Matched  Filter 
on  Raw 
Seismogram# 

S/N 

Matched  Filter 
on  Band  Paaa 
Filtered 
Seismograms 

B1 

3.1 

3.4 

8.3 

8.6 

P3 

i.e 

4.7 

10.2 

10.3 

n 

3.3 

4.5 

8.5 

8.7 

AO 

3.6 

3.9 

9.0 

9.0 

B3 

4.3 

4.4 

9.1 

9.0 

C4 

3.4 

4.1 

9.2 

9.4 

B4 

3.2 

3.8 

7.9 

7.8 

Cl 

4.0 

3.5 

7.7 

7.6 

C2 

4.1 

4.4 

9.7 

9.7 

B2  '* 

4.6 

4.5 

9.1 

9.1 

C3 

4.7 

4.7 

8.7 

8.7 

D3 

3.7 

4.0 

9,0 

9.0 

D4 

4.3 

4.6 

7.0 

7.1 

Dl 

3.7 

4.4 

10.5 

10.3 

02 

4.9 

6.0 

11.0 

11.0 

*3 

4.3 

4.6 

8.4 

8.5 

E4 

4.6 

4.4 

9,4 

9.4 

El 

3.7 

3.3 

4.7 

4.6 

FI 

3.3 

3.3 

6.4 

7.0 

E2 

2.6 

2.6 

4.9 

5.0 

F2 

3.3 

2.9 

5.4 

5.7 

Table  III-B . 

Observed  S/fo  ValueH  for  LRSM  Stations 


STATION 

B/N 

Raw 

Beiaaograma 

B/N 

Band  Pfisa 

filtered 

Seiamograma 

8/N 

Hatched  Filter 
on  Raw 
Selawograaa 

B/N 

Hatched  rilter 
on  Band  Peaa 

Fi lterad 
Beiamograma 

SV3QB 

2.4* 

3.2* 

1.9 

2.i 

HN-ME 

3.0* 

3.4* 

1.8 

2.3 

HK-ON 

7.5 

9.6 

14.3 

14.7 

po-ac 

B.  9 

16.3 

17.1 

17.3 

KC-MO 

6.4 

6.0 

9.2 

8.8 

BHO 

4.4* 

3.7 

7.8 

9.3 

MO-ID 

C.l 

4.9 

10.5 

10.2 

UBO 

6.5* 

3.9 

4.5 

7.7 

AX2AL 

3.1 

6.7 

10.0 

10.4 

EU2AL 

1.8 

3.9 

5.7 

6.2 

HMO 

6.3 

8.6 

15.4 

15.2 

M-rt 

3.4* 

3.5* 

2.2 

2.6 

JC-LA 

2.0* 

4.6* 

1.0 

4.5 

MN-HV 

3.3* 

3.8* 

2.2 

2.2 

*  Hole*  Peak  Hi ■ taken  For  Signal  by  computer 


Predicted  S/S  Gs-in 


Predicted  S/N  Gain 


RELATIVE  MAGNIFICATION 
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Figure  2.  LASA  Long  Period  System  "A"  Response 
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Figure  4.  LRSM  Long  Period  System  Response 
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Figure  5.  Digital  Band-Pass  Filter  Response 


Figure  6.  Data  Processing  Flow  Chart 


Several  Signal  Enhancement 
■  A  Single  LASA  Element 


Figure  8.  Comparison  of  Several  Multichannel  Signal  Enhano 

Procedures  for  21  LASA  LPZ  Elements 


i-Channel  Combinations  of  9  LASA  LPZ  Elements 
Effect  of  Inter-Sensor  Spacing 


11.  Typical  Results  of  Matched  Filter  Processing 
of  LRSM  (s  VELA  Observatory)  LPZ  Seismograms, 
and  the  Sum  of  14  Such  Matched  Filter  Outputs 
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Figure  14.  Rayleigh  Wave  Enhancement  vs.  Inter-Sensor  Spacing  (N=9) 
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Both  singls  chsnnel  and  array  signal  enhancement  tech- 
nlquee  hav*  been  applied  to  Rayleigh  waves  from  a  small  Green- 
land  Saa  earthquake  recorded  at  LA8A  and  13  LRSM  or  Observatory 
stations.  Analysis  of  individual  LASA  lonq  period  recordings 
indicated  that  the  matched  filter  increased  S/N  by  more  than 
*  db  over  the  mean  S/N  of  band  pan  filtered  (15-50  sec  period) 
seismograms.  Band  paae  filtering  increased  the  mean  S/N  by 
only  1.5  db  over  the  mean  for  unfiltered  seismogramo  for  LASA. 

Additional  improvement  in  8/N  from  beam-forming  both  band 
passed  and  matched  filtered  traces  approached  the  expected  in¬ 
crease  for  uncorrelated  noise  provided  the  intersenaor  spacing 
(mesh-siie)  was  at  least  30  km.  Comparable  improvement  was 
obtained  for  both  single  channel  and  beam- formed  LRSM  data, 
for  LASA  tha  beam-formed  matched  filter  8/N  for  13  stations 
was  17  db* above  the  mean  8/N  for  tha  individual  band  pass 
filtered  aeiamogramaf  that  for  13  LRSM  stations  was  151*  db 
sbova  the  mean  S/M  of  band  paes  filtered  seismograms.  Beam¬ 
forming  matched  filter  osismograsia  consistently  produced  s/N 
vsluss  7-9  db  above  the  s/N  for  beam -tf-ed  band  pace  filtered 
emismogrene.  The  effects  of  such  arre<>  parameters  as  number 
of  sensors#  sensor  spicing#  and  aperture  on  signal  enhancement 
ere  evaluated  for  this  event. 
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